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Abstract. It has long been accepted that a possible mechanism for explaining the 
existence of magnetic white dwarfs is the merger of a binary white dwarf system, as 
there are viable mechanisms for producing sustainable magnetism within the merger 
product. However, the lack of rapid rotators in the magnetic white dwarf population 
has been always considered a problematic issue of this scenario. In order to explain 
this discrepancy we build a model in which the interaction between the magnetosphere 
of the star and the disk induces angular momentum transfer. Our model predicts that 
the magnetospheric interaction of magnetic white dwarfs with their disks results in a 
significant spin down, and we show that the observed rotation period of RE J 03 17-853, 
which is suggested to be a product of a double degenerate merger, can be reproduced. 



1. Introduction 

Magn etic white dw arfs (MWDs) — white dwarfs with field strengths rang ing from 
1 kG (IJordan et al.ll2007h up to 1 GG dKawka et al.ll2007l : iKiilebi et al. 2009J) — com- 



prise more than -10% of all white dwarfs. There are two possibilities to account for the 
observed magnetic fields. According to the fossil field hypothesis, these white dwarfs 
descend from Ap/Bp stars ( Angel et ai]|l98lh . and their magnet ic fields are remnants of 



the previous evolution. However, population synthesis studies (IWickramasinghe & Ferrario 



2005) show that the observed number of Ap/Bp stars is insufficient to explain the 
number of MWDs, as well as other properties of their distribution. One mechanism 
which has been suggested as viable for producing ultra-massive white dwarfs (Mwd > 
1 . 1 Mq) is the merger of two white dwarf s, which naturally explains why MW Ds are 



more mas sive than their field c ounterparts (Mars het al.ll 19971 : iKepler et al.ll2007 ). 



iGarcfa-Berro et ajj (120121) proposed that these magnetic fields are generated in the 
corona above the merger product through an aa> dynamo, since convection and differen- 



tial ro tation exist simultaneously. Using the results of the simulations of lLoren-Aguilar et al 



(2009) they showed that very powerful magnetic fields can be generated in the hot 



corona resulting from a merger. Also, to account for the slow rotational velocities 



1 



2 



B. Kiilebi et al. 



of the MWDs, they suggested magneto-dipole radiation as a possible solution for this 
drawback of the model. However, it is likely that magneto-dipole radiation would not 
be strong enough to slow down the remnants of the coalescences to the observed rota- 
tion periods of MWDs. Nevertheless, if the merger product has high magnetic dipole 
fields and the disk from the merger survives, it is possible that the long-term evolution 
of the system will involve angular momentum exchange between the disk and the star, 
depending on the magnetic field strength of the central object. In this work we show 
that this angular momentum exchange successfully spins down the central MWD to the 
observed periods. 



2. Debris disk 



To model the accretion disks which are remnants of double degenerate mergers, we con- 
sider a freely expanding thin disk model which is similar to the disks fo rmed by tidal 
disruption of stars by super-massive black holes ( Cannizzo et al.lll99Ch or supernova 
fallback disk s around young neutron stars, e.g. anomalous X-ray pulsars dChatterjee et al 
120001 : lAlparlfeOOll) . In this type of model, the ev olution of the disk can be described in 
terms of its initial mass and angular momentum (lErtan et alJl2009h . 

The temporal evolution of t he surface mass density in the disk is described by a 
diffusion equation (Pringlel ll98lh : 
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where £ is the surface mass density an d v is the turbule nt kinematic viscosity. The 
rest of the thin disk st ructure equations (|Frank et al.ll2002r) can be solved algebraically 
( Cannizzo et al.|[l99fj|) to obtain a viscosity in the form v — Cr p lL q where C , p and q 
are constants determined by the opacity regime — see also lErtan et all (12009b . For this 
type of viscosity, Eq. (Q} has the self-similar solutions found by lPringlel (11974) . One of 
these solutions conserves the angular momentum of the disk (/d = 0), and the mass of 
the disk evolves as a power-law in time. In particular, for this solution the mass flow 
rate at the inner boundary can be expressed as: 



M = M 1 1 + - 

to 



M = (a- 

to 



(2) 



where to is the viscous timescale at some characteristic radius, Mo the initial mass of 
the disk and Mq is the initial accretion rate. The power-law index a c an be written in 



terms of p and q, which in turn are determined by the opacity regime (|Cannizzo et al. 
1990). For electron scattering the exponent is a = 19/16 and for the bound-free opacity 
regime is a = 5 /4. 

The most important parameter governing the spin evolution of MWDs is the vis- 
cous t imescale to, whic h is determined by the initial mass and angular momentum of the 
disk (lErtan et alJliOfflh . For our calculations we adopt a metal disk in the bound-free 
opacity regime with a S s = 0.1, for which the viscous timescale is 



MWDs with debris disks 



3 



where jo = Jq/Mq is the specific angular momentum. Note that the viscous timescale 
depends strongly on jo- 

The way in which the star interacts with the disk depends on the interplay between 
the angular velocity of the star Qwd and the keplerian velocity at the inner radius of 
the disk Ok- This is parametrized by the fastness parameter w» = Qwd/^kWiiX from 
which the torque transfer can be computed: 

N d = n(co*)N . (4) 

where A^n = frc CRin ) R \ M and «(a»*) is called the dimensionless torq ue dGhosh & Lamb! 
19791 : lLamblll988L (l989h . Following lAlparl (l200lh and lEksi et all (|2005|) we employ 
n - 1 - a> t . If the inner radius of the disk goes beyond the light cylinder radius, the disk 
can no longer torque the star via the magnetosphere, = 0. In this so-called ejector 
stage the star will spin-down via magneto-dipole radiation 



3. Results 

To compare our model with observations, we first conducted a set of simulations in 
which the adopted masses of the central remnant and the disk, the angular velocity of 
the uniformly rotati ng star and the angular mo mentum of the disk were those of the 
SPH simulations of lLoren-Aguilar et all (12009b . For a given simulation, we calculated 



the spin evolution of the MWD for different polar field strengt hs. To determine the 
total integration time we used the simple cooling law of iMesteJ (U965) and estimated 
the cooling age at the otherwise typical observed value of T e g~ 40000 K. The exact 
value of of the initial core temperature depends on the mass of the merging components. 
However, the exact value of the initial core temperature has virtually no consequences 
for the age estimate, as white dwarfs cool very rapidly during the early phases of evo- 
lution. Under these conditions the cooling ages for MWDs of masses 0.8, 1.2, 1.4 Mq 
MWDs were 13, 50 and 100 Myr respectively. 

The evolution of the angular velocity of the star is determined by the fastness 
parameter, a»*. Initially the accretion rate is very large. Hence, the magnetospheric 
radius is very small, and the inner radius of the disk is equal to the stellar radius. In 
this case, the fastness parameter is smaller than 1. Thus, the remnants of mergers are 
fast rotators, and are rotating at almost one third of the critical angular velocity of the 
star co t ~ 0.33 s . As time passes by, angular momentum is transferred onto the star 
and the fastness parameter approaches 1. If the mass of the white dwarf is not close 
to Mch the star does not reach critical rotation and this spin-up stage continues till 
<y* > 1. At this point the magnetic torque from the disk starts to spin down the star. As 
the accretion rate decreases, the magnetospheric radius and hence a;*, increase rapidly. 
This results in a stronger torque on the star which eventually f orces it to reach torqu e 
equilibrium at a** - 1. This phase is called the tracking phase dChatterjee et alj |2000). 
This is a quasi-equilibrium phase because the mass transfer rate keeps declining and the 
equilibrium period, at which <y* ^ 1, increases. Hence, accretion continues while the 
star spins down. This second phase of torque transfer is observed in all our simulations 
and the strength of the magnetic field determines when this stage is reached. If the 
polar magnetic field B p is strong, the Alfven radius and oj* are also larger. In all our 
simulations there is a limiting magnetic field strength which determines whether an 
MWD will enter this strong spin down phase or not, given the same cooling age. If 
the magnetic field is close to this limiting value, Bi; m , the final period of the star is 
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Figure 1. The evolution of the period for a 1.30 Mq MWD for a turbulence turn- 
off temperature T p = 1000 K. The dashed lines show where the torque transfer from 
the disk becomes ineffective and magneto-dipole spin down takes over. 



very sensitive to the exact value of B p . In particular, for a field strength of B \\ m , a 
10% difference in the initial field strength results in a change by a factor ~ 2 of the 
final rotation period. Finally, after this strong spin down phase, if the magnetic field 
strengths and cooling age are adequate, the system reaches quasi-equilibrium around 

co* < 1- 

The period reached after the tracking phase depends on time, dipole field strength, 
angular momentum of the disk, mass and radius of the star, but not on the mass of the 
disk. Building upon this, and taking into account that there are considerable unknowns 
about the super-Eddington accretion phase, we also built models with low-mass disks. 
This approximation is acceptable since this phase of the evolution is short (10 4 yr) with 
respect to the evolutionary time of the observed population of MWDs (10 7 - 10 8 yr). 
For these models we assumed a disk mass of 1(T 3 M Q — the mass of the disk when 
the Eddington accretion phase finishes and accretion driven by the evolution of angu- 
lar momentum takes over — while the rest of the material was incorporated onto the 
remnant white dwarf. We performed this test for a 1.2 Mq MWD, and the results were 
identical to those of obtained in the previous set of simulations, for cooling ages longer 
than ~ 10 7 yr and B p > Si im ~ 10 MG, which apply for most high-field MWDs with 
known periods. 

The lifetime of a MWD is expected to be long when compared with the typical 
cooling timescale of the disk. Thus, it is natural to wonder whether or not turbulence 
in the disk can be sustained during long periods of time. The cooling t imesc ale of the 



disk depends primarily on its degree of the ionization. iMenou et al.l (120011) showed 
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that for thin disks made of C, thermal equilibrium is reached near 5 000K, and that 
for all compositions recombination always occurs below ~ 1 000 K. Although this tem- 
perature is high, it is known that these disks can be partially ionized at even lower 
temperatures due to severa l other processes , like cosmic ray heating, collisions and 
charged-particle mobility. Ilnutsuka & Sand (l2005h carefully studied all these effects 
and reached the conclusion that at temperatures < 1 000 K, magneto-rotational instabil- 
ity dBalbus & Hawlevlll99lT) could be still effective. 

To implement the effects of the turbulent life time, we calculated the temperature 
at the inner radius of the disk at each time step. If this temperature was smaller than 
T p - 1 000 K, we assumed that the disk was completely inactive and accretion was 
switched-off. Thus, the MWD immediately enters in the so-called ejector phase and 
spins down only by magneto-dipole radiation. To compare these calculations to the case 
of RE J 0317-853 — wi th a period of 721 s dBarstow et alJI 1995b iFerrario et alJll997h 
— and WD 1658+441 dSchmidt et all 1 19921) we applied this procedure to the model 
with M = 1.30 Mq. The evolution of the periods for MWDs with different dipole 
fields is shown in Fig.Q] We found that although the magneto-rotational instability is 
long-lived (~ 10 5 - 10 7 years), the disk becomes inactive at the current ages of these 
objects (~ 10 8 years). Moreover, since the magneto-dipole torque is not as strong as the 
magnetic torque from the disk, the final period is determined by the time at which the 
magneto-rotational instability is turned off, which depends on the accretion rate. 



4. Conclusions 

For MWDs with field strengths between 50 MG and 1 GG, and T p ranging from ~ 
300 K to 1000 K, the final periods reached are around 10 3 - 10 4 s (see Fig.Q]), in 
agreement with the period of RE J 0317-853. However this is not true for the MWD 
WD 1658+441 which has a field s trength of « 2MG and undetected — hence, a long 
(> 10 5 s) — period dShtol' et alJll997h . RE J 0317-853 has been suggested to be a 



merger remnant due to its mass, cooling age and strong magnetism. Our simple evo- 
lutionary models of the coupling between the white dwarf and the surrounding disk 
show that the disk may remove angular momentum from the central white dwarf for 
reasonable choices of the free parameters, and thus these models are able to explain the 
observed properties of REJ 0317-853. 
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